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O '■ Abstract 
in 

' The apparent excess of Z — > 66 events at LEP may be an indication of new 

a: 

O 1 physics beyond the standard model. However, in either the two-Higgs-doublet 

■ model or the minimal supersymmetric standard model, any explanation would 



lead to an important new decay mode of the top quark and suppresses the 



5— 1 

■ t — ► branching fraction, which goes against what has been observed at the 

Tevatron. In the two-Higgs-doublet model, the branching fraction of Z —> bb 
+ a light boson which decays predominantly into bb would be at least of order 
10" 4 . 
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I. INTRODUCTION 



With an accumulation of 8 x 10 6 Z decays into hadrons and charged leptons by the four 
LEP experiments at CERN by the end of 1993, measurements of a large number of rates, 
branching fractions, and asymmetries have now become even more precise 0. The only 
apparent deviation by two or more standard deviations from the prediction of the standard 
model is in the ratio 

R b = T(Z - hl) (1) 
T(Z — > hadrons) 

Assuming mt = 175 GeV and m# = 300 GeV, the standard model predicts that Rb = 0.2158, 
whereas LEP obtained Rb = 0.2202 ± 0.0020 if the similarly defined R c is assumed to be 
independent. If the latter is fixed at its standard-model value, then Rb = 0.2192 ± 0.0018. 
In either case, the excess is about 2% ± 1%. If this is taken seriously, physics beyond the 
standard model is indicated. 

In this paper we will examine two frequently studied extensions of the standard model: 
the two-Higgs-doublet model (2HDM) and the minimal supersymmetric standard model 
(MSSM). We will assume that the only significant deviation from the standard model is Rb, 
hence we will take Rb = 0.2192 ± 0.0018 as the experimental value and see how these two 
extensions may be able to explain it. We will concentrate on obtaining Rb > 0.2174, i.e. 
within one standard deviation of the experimental value, because otherwise the difference 
with the standard model is insignificant and we might as well not bother with any possible 
extension. 

Whereas the contributions to Rb from either the 2HDM |2|-[5| or the MSSM have 
been studied previously, we are concerned here also with the effect of these new contributions 
on top decay. Together with the constraints from the oblique parameters |J, we find that 
a large Rb excess will always lead to an important new decay mode of the top quark and 
suppresses the t — > Wb branching fraction, which goes against what has been observed at 
the Fermilab Tevatron [|10||. In the 2HDM, we also find that the branching fraction of Z — > bb 



+ a light boson which decays predominantly into bb will be at least of order 10 pi 



II. TWO HIGGS DOUBLETS 



The simplest extension of the standard model is to have two Higgs doublets instead of 
just one. The relevance of this model to Rb was studied in detail already a few years ago 
0. To establish notation, let the two Higgs doublets be given by 

\ 













\ 



<f>: 



+ 



K (vi + rii + iXi)/V2 j 



(2) 



Let tan/5 = v 2 /v 1 , then 



h + = <p2 cos (3 — <f>f sin j3 
A = Xi cos/3 - xi sin/3 , 
hi = rji cos a + r\ 2 sin a , 
h 2 = T)2 cos a — r)i sin a . 



(3) 
(4) 
(5) 
(6) 



The corrections to the left- and right-handed Zbb vertex induced by the charged Higgs 
boson h + and the neutral Higgs bosons hi, h 2 and A are given by B 



a 



Ans 2 w 



\2 C I 2 2 2\ 

K,r Pl,r{1 > m +> m J 



(7) 



and 

6gl R (hi^A) 
where 



a 



' m b tan/3" 



r 2 Nil 2 2 \ , • 2 Nil I 2\l /o\ 

[cos a p LR (g ,m l5 m A ) + sin a p LijR (g ,m 2 ,m A )J , (8) 



p£ o(g 2 , m 2 , M 2 ) = {-, +}> 2 C (g 2 , m 2 , M 2 , M 2 ) + g* p 3 {q\ m 2 , M 2 , M 2 ) 



p 3 (q 2 , m\ M , M ) - p^q\ m\ m\ M ) 



pl R {q\ m\ M ) = {-, +}p A (q\ m\ M , 0) + g> p 3 {q\ m\ 0, 0) + p 3 {q\ M , 0, 0) 



(9) 
(10) 



In the above, g L = — \ + \s"w, g R = gS 2 ^ and g h = \ — s^r. We have also assumed that 
$2 couples to the up-type quarks and $i to the down-type quarks, hence = cot (3, 



and Xr = y ^ lw tan p. The masses of h + , hi, h 2 and A are denoted by m + , m%, m 2 , and 
rriA respectively. The functions C and p 3j4 are defined in Ref. [2]; see also Appendix 

For a heavy top quark, it is well-known that pi R (q 2 ,m 2 + ,m 2 ) ~ {— ,+}(— |). Hence 
g b L Re{6g b L (h + )} and g b R Re{6g R (h + )} are negative, thereby decreasing the value of Rb. This 
means that the tan (3 < 1 region can be ruled out In this region, t — > bh + also 

becomes the dominant decay for the top quark |l2j unless it is kinematically not allowed. 
In fact, any significant reduction of the t — > Wb branching fraction is in conflict with the 



Tevatron data [jlOfl because the number of top events observed is such that even if we assume 
B(t —>■ Wb) = 1, the deduced experimental ti production cross section is already larger than 
expected |13| . 

If tan/3 is large, the contribution from the neutral Higgs bosons becomes important. In 
other words, Eq. (||) must be considered even though it is suppressed by (mb/m w ) 2 . Note 
that since Rb is proportional to {g b L + 5g b L ) 2 + (g R + Sg b R ) 2 , p^ is more important than p R 
because g\ » g R 0. Again because g R is small, p^ is dominated by p 4 (g 2 , m 2 , M 2 , 0) in 
Eq. (10). In order that p^{q 2 , m 2 , M 2 , 0) be positive and not too small, both m and M must 
be light, namely m 2 ,M 2 < q 2 . In particular, for Pi{m 2 z , m 2 , M 2 , 0) > 0.2, both m and M 
should be less than 65 GeV. 

It was shown already in Ref. [2] that for tan/3 = 70 ~ 2m t /mb, the Rb excess peaks at 
about 4% near = mi ~ 40 GeV for a = 0. However, since Z — > Ahi is not observed, 
mA + mi > mz is a necessary constraint. We show in Fig. [I] the contours in the mi — mA 
plane for Rb = 0.2192 and 0.2174. It is clear that relatively light scalar bosons are required 
if the Rb excess is to be explained. 

For A(hi) lighter than mz and having an enhanced coupling to bb, the decay Z — > 
bb + A(hi) becomes nonnegligible [13]. As an illustration, we show in Fig. |^ the branching 



fractions of these two decays as functions of with the constraint + mi = mz + 10 
GeV so that a reasonable fit to the Rb excess is obtained. It is seen that the sum of these two 
branching fractions is at least of order 10~ 4 . Once produced, A or hi decays predominantly 
into bb as well. Hence this scenario for explaining Rb can be tested at LEP if the sensitivity 



for identifying one bb pair as coming from A or hi in bbbb final states can be pushed down 
below 10" 4 . 

Since &£ is involved in any enhanced coupling to light particles in explaining the Rb excess, 
its doublet partner t^ must necessarily have the same enhanced coupling to related particles. 
In the 2HDM, we must have an enhanced ibh + coupling. Therefore, unless m + > m t — m^, 
the branching fraction of t — > bh + will be important. As a result, the standard t — > Wb 
branching fraction will be seriously degraded. We show this in Fig. |as a function of m+. 
Large values of m + are disfavored in this scenario because the splitting with A and hi 
would result in a large contribution to the oblique parameter T, resulting in the constraint 
m+ < 150 GeV g. 



III. SUPERSYMMETRIC HIGGS SECTOR 

In the minimal supersymmetric standard model (MSSM), the two Higgs doublets have 
exactly the same gauge and Yukawa couplings as in the 2HDM we discussed in the previous 
section. In addition, the quartic scalar couplings of the MSSM are determined by the gauge 
couplings and there are only three arbitrary mass terms: $2 $2, and <&\<&2 + ^2^1 • 

Hence we need only two extra parameters, usually taken to be tan/3 and m^, to specify the 



entire Higgs sector, subject of course to radiative corrections [15|. However, these corrections 
are only significant for small tan (3 and since we need a large tan /3 to explain Rb, we will use 
the simpler tree-level expressions in our numerical analysis. Combining both the charged 
and neutral Higgs contributions, we plot in Fig. [| the Rb = 0.2192 and 0.2174 contours in 
the rriA — tan/3 plane, with the constraint + mi > mz- We plot also the contours for 
B{t — > Wb) = 0.85 and 0.7, which correspond to reductions of 28% and 51% of the top 
signals at Fermilab respectively. It is abundantly clear that the Higgs sector of the MSSM 
is not compatible with both a large Rb and a small B{t —>■ bh + ). In the 2HDM, tjia and 
m + are independent parameters, whereas in the MSSM, there is the well-known sum rule 
m 2 + = m 2 A + mw- Hence an approximate custodial symmetry exists in the MSSM to keep the 
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contribution to T small, but at the same time m + < m t — m;, is inevitable if is assumed 
to be small enough to obtain a large Rb excess. In Fig. |5| we plot the minimum branching 
fraction B(Z — > bbA + bbhi) for a given lower bound of Rb- This minimum is obtained by 
varying tan/? with a fixed value of tua- As in Fig. 2, we see that this branching fraction is 
at least of order 10~ 4 . 



IV. CHARGINOS AND NEUTRALINOS 

In addition to the Higgs contributions, there are chargino (x) and neutralino (M) con- 
tributions to Rb in the MSSM. They have been studied previously ||-|§, and it is known 
that the particles in the loops have to be light in order to obtain large contributions. The 
parameters involved here are tan/3, //, M 2 , ^ 12 , m b 12 i an d 9 b , where the supergravity 
condition Mi ~ O.5M2 at the electroweak scale has been assumed. The scalar mixing angles 
Q l and 6 h as well as others are defined in Appendix [B]. 

Since scalar quarks have not been observed at LEP, their masses must be greater than 
half of the center-of-mass energy, namely g > \mz- On the other hand, the lightest 
neutralino (Mi) is always lighter than the lightest chargino (xi)> thus the condition m xi > 
\mz is not enough by itself. Let us define the conservative constraints from the invisible 
width and total width of Z as 5T mv = r inv (expt)| 

max (SM)| min and ST Z = r z (expt)| max 
- r^(SM)| min , where SM denotes the standard model. 
From the updated LEP data we obtain 

T(Z -> M1M1) < 5T inv = 7.6 MeV , (11) 
T(Z -»• MiMj) < 8T Z = 23 MeV . (12) 

These constraints must be included for the analysis in order to provide consistent results. 

The corrections to the left- and right-handed Zbb vertex induced by the charginos (x) 
and neutralinos (M) are given by 

5gi R (x,M) = ^-F L>R (x,M) , (13) 
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where 



-\js, A iA .' p 4 (g ,m t -.,m ?j .,m x J - A fci ' Cy A fc j p 3 (q ,m h ,m x .,m x .) 
+ A ki R (£>5 ,L - Oij R ) A kf* m xr m xj C (q 2 , m| fe , mj. , m^. ) , (14) 
F l ,rW = X jk R * S ji X ik R P^1 2 ^ m t^ m t^ m M k ) ~ X kf Qij R X kf* P3(q 2 ,ml k ,m 2 K ,m 2 ^.) 

+\l {Qij - Qy' ) X k j m Ni m Nj C {q , m-^ m H ., m N .) , (15) 

and there is an implicit sum over all repeated indices. See Appendix B for the definitions of 
the various quantities in the above. 

For tan/5 < 20 or mg > mz, the chargino contribution is the most important. From 
Eq. (|I~4"D , we see that the contribution is the largest when the lighter scalar quark t\ is mostly 
tji, namely Q\y = 0. We follow the usual strategy of finding the maximally allowed for 
a given tan f3 and m Xl . Here we impose also the LEP constraints given by Eq. (|I~2"D . In Fig. |6|, 
we plot the maximally allowed Rb as a function of tan/3 for m Xl = 60 GeV and = \mz- 
We see that R^ > 0.217 can be obtained. However, top decay into t\ and a neutralino is 
now possible and the corresponding branching fraction B(t —> Wb) shows clearly that this 
solution would conflict with the Tevatron data [ PH| , P^ |. 

If tan/3 is large, the neutralino contributions become important because the b quark 
coupling to the higgsino is proportional to 1/ cos/?. Here we would like to point out that our 



biMf b couplings given in Eqs. (|B15|) and (|B16j ) are different from those given in Ref 



but agree with Ref. |16]. For simplicity and with little loss of generality, we assume that 
m~ tx = m~ bi = 60 GeV and m^ 2 = m~ b2 = 250 GeV. We vary the scalar quark mixing angles 
9 l and 8 b , such that Rb is maximum within the allowed LEP constraints given by Eq.(|T2"D 
and m xi > \mz- This is the most optimistic scenario; we cannot achieve a large enough Rb 
otherwise. Taking tan/? = 70, the contours of Rb = 0.2174 and 0.2192 are plotted in Fig. 
0. Again we plot the t — > Wb branching fraction. We find only very narrow regions where 
B(t — > Wb) > 0.7 and Rb > 0.2174. Hence future experiments on top decay will play a 
decisive role to verify or rule out this scenario JT7]]. 

The dominant contribution to the oblique parameter T comes from the scalar quarks. 



We have checked that T ~ 0.4 in the narrow regions, which is in mild conflict with the 
recent global fit T = — 0.67± 0.92 [IB]. Nevertheless, Fig. [7] represents the most optimistic 
scenario. In Fig. |8j, we consider a more restrictive case with 9 l = 9 b so that T = 0. As a 
result, the narrow regions shrink as expected. 



V. CONCLUSION 

The 2HDM p|-^J and the MSSM have each been suggested to explain the Rb excess 
at LEP. Here we consider also the effects of these new contributions on top quark decay. In 
the 2HDM, large tan/3 and light neutral scalars are necessary to increase Rb to within one 
standard deviation of the experimental value. The corresponding ibh + coupling then allows 
top decay into b and h + unless it is kinematically not allowed. However, m + > 150 GeV 
would be in conflict with the constraint of the oblique parameter T. The same interactions 
which allow a large Rb also allow the decays Z — > bbA{hi). We show in Fig. 2 that the 
branching fraction B(Z — > bbA + bbh\) is at least of order 10 -4 which can be tested at future 
LEP experiments. 

In the supersymmetric Higgs sector, there are only two independent unknown parameters: 
tan/3 and m^. Specifically, because of the sum rule m+ = m\ + m^, top decay is always 
possible for a small enough rriA to account for the Rb excess. We see in Fig. f| that the 
region which allows Rb to be large does indeed conflict with top decay. In addition, B(Z —> 
bbA + bbhi) is also at least of order 10~ 4 as shown in Fig. |5|. 

Because of the scalar-quark mixing angles, the chargino and neutralino contributions to 
Rb could each be either positive or negative. Here we consider the most optimistic scenario 
that the mixing angles are chosen to maximize the Rb value. The chargino contribution, as 
opposed to the charged Higgs and W contributions, can increase Rb above 0.2174 if t\ and 
Xi are light enough. Since M\ is always lighter than xi, we impose also the LEP constraints 
given by Eq. flT2p. Again this new contribution gives rise to a new channel for top decay 
which reduces the t — > Wb branching fraction significantly. For large tan /3, we consider both 



S 



chargino and neutralino contributions. We find that there are only very narrow regions, as 
shown in Figs. [7] and 8, in which both large Rb and B(t — > Wb) are compatible. Future 
experiments on top decay would verify or rule out this scenario. 

Note added. After the completion of this paper, additional data from LEP have been 
reported. [I9~| For R c fixed at its standard-model value, the latest experimental value of Rb 
is 0.2205 ± 0.0016. This means that our assertion that Rb conflicts with top decay in the 
2HDM and the MSSM becomes even stronger. 
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APPENDIX A: DEFINITION OF p 3 



The general form for p 3 is given by 



p 3 (g 2 ,m 2 ,M 1 i ,M 2 2 ) = [q 2 (C 22 ~C 23 )+2C 24 -M 1 M 2 C ](q 2 ,m 2 ,Mf,M 2 ) 

+ ^[B 1 (0,M 2 ,m 2 )+B 1 (0,M 2 ,m 2 )] - i. 



(Al) 



For Mi = M 2 , P3 is identical to that defined in Ref. 0. This generalized definition is useful 
in the supersymmetric case. 

APPENDIX B: MASSES AND MIXING IN THE MSSM 



In this appendix, we present the relevant couplings for the chargino and neutralino 
contributions to the Zbb vertex. We allow general scalar top mixing, namely ti = 0\i ti 
and tR = ti with a similar definition for the scalar bottom quarks. The Zt*tj and Zb*bj 
vertices are given by S\f(pi - pj)^, where 



4 = ^X-^in 2 W„ 
I „i * „i 1 



WOij 



S 



-9 b u ^. + ±sin 2 ^ 



(Bl) 
(B2) 



The p's are defined as out-going momenta. 

The chargino and neutralino mass matrices are given by 

M 2 \plm-w cos (3 

\[2m w sin (3 [i 

which links (iW~, h{) T to (iW + , h%), and 
/ Mi 

M 2 
— mz sin 9 W cos (3 m z cos 6w cos /3 
\ sin 6*vi/ sin j3 —mz cos 6*^ sin (3 



(B3) 



— sin 9w cos /? 
m z cos 9w cos /? 




m z sin sin /? \ 
—m z cos 6V sin (3 
-H 







(B4) 
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in the (iB ,iW s , h®, h®) basis. The mass eigenstates, x% an d -M, are related to these basis 
states by the following transformations: 

(iW + , ht f = V ijXj ; (iW~, K) T = Uvxf , (B5) 

and 

(iB, iW 3 , til, ti° 2 ) T = NijUj . (B6) 
Thus, M x and are diagonalized by 

V T M x U = m Xi 8 i3 , (B7) 

and 

N T M M N = mjsrAj ■ (B8) 

Let us first consider the vertices for loops involving charginos x- The ZxJXj vertex is 
g^en by ^7,[^ ^ + Og ^] with 

Q% = cos 2 6 W 5 ZJ - ±V 2 * V 2j , (B9) 
Ofj = cos 2 9 W Sij - \U 2i U% , (BIO) 

and the fj vertex is given by g[Afj + Ag ^p 2 -] with 

A S = ^ + ^ - - V* ^ , (Bll) 

Ag = [/* 9{* . (B12) 

y2m w cos p 

For the vertices involving neutralino loops, ZTT-Afj vertex is given by cos g g w 7/JQg ^p- + 



Qg ^] with 



qS = I^^-^^)» ( B13 ) 

Qg = \(N 4l N* i3 - N M 7V|.) , (B14) 



and the 6; vertex is given by g[Ag ^p- + Ag with 
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A 



V2 



— tan 6% J\ u V 2j -= 



\f2mw cos j3 

i * 

* no 



N* 6 



3i u lj 



(B15) 
(B16) 



iw cos /3 

Note that the relative signs of the m terms in Eqs. ( P15| ) and (|B16| ) are in agreement with 



Ref. TO, but differ from Ref. 
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FIGURES 

FIG. 1. Rb = 0.2192 (solid) and 0.2174 (dashed) contours in the m\ — tua plane for a = 
and tan/3 = 70. The straight line corresponds to tua + mi = Mz- We have also assumed 
m + = ?7i2 = 175 GeV. 

FIG. 2. The branching fractions, B(Z -> bbA) (dashed) and B(Z -»■ 66/ii) (dotted) and their 
sum (solid), as functions of wia where we take via + m\ = Mz + 10 GeV, tan (3 = 70, a = 0, and 
m + = 777,2 = 175 GeV. 

FIG. 3. The branching fraction B(t — > W6) as a function of m+ for tan/3 = 70 (solid), 50 
(dashed), and 20 (dotted). 

FIG. 4. i? 6 = 0.2192 (heavy), 0.2174 (solid) and B(t -► W6) = 0.85 (dashed), 0.7 (dotted) 
contours in the 777,^4 — tan/3 plane. We also plot the constraint 777,4 + mi > mz (dash-dotted). 

FIG. 5. The minimum branching fraction B{Z — > + 66/ti) for a given lower bound of R^, 
where we take mA = 45 GeV (solid) and = 50 GeV (dashed). 

FIG. 6. The maximally allowed Rf, (solid) as a function of tan (3 for m xi = 60 GeV and 
m ti = \ m z- We also plot the corresponding branching ratio B(t — > bW) (dashed). We have 
assumed 777^ = 250 GeV and 0\ 2 = 1. 

FIG. 7. Contours of maximally allowed values Rf, = 0.2174 (solid) and 0.2192 (dotted) as well 
as B(t — > IF6) > 0.7 (dashed) in the \i — M2 plane where the heavy lines represent the LEP 
constraints and 777 Xl > \mz- We assume 777^ = = 60 GeV and = m^ = 250 GeV. 

FIG. 8. The same as Fig. 7 but with the added constraint Q l = b so that T = 0. 
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